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ABSTRACT 

We present a unified derivation of the resummation of Sudakov logarithms, directly from 
the factorization properties of cross sections in which they occur. We rederive in this manner 
the well-known exponentiation of leading and nonleading logarithmic enhancements near the 
edge of phase space for cross sections such as deeply inelastic scattering, which are induced 
by an electroweak hard scattering. For QCD hard-scattering processes, such as heavy-quark 
production, we show that the resummation of nonleading logarithms requires in general mixing 
in the space of the color tensors of the hard scattering. The exponentiation of Sudakov 
logarithms implies that many weighted cross sections obey particular evolution equations 
in momentum transfer, which streamline the computation of their Sudakov exponents. We 
illustrate this method with the resummation of soft-gluon enhancements of the inclusive Drell- 
Yan cross section, in both DIS and MS factorization schemes. 
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1. Introduction 



Our aim in this paper is to treat cross sections and amplitudes that allow Sudakov ex- 
ponentiation and resummation in a unified and streamlined manner. We shall point out the 
factorization that underlies these resummations, and derive the resummed form directly from 
it. The close relation of factorization to Sudakov exponentiation was emphasized and ex- 
ploited long ago by Mueller) 11 Collins and Soper' 21 and Sen! 31 Here, we shall generalize this 
observation, and develop a method that makes the wide range of its application to cross 
sections and amplitudes straightforward and, we hope, clear. 

A classic example W of Sudakov exponentiation is the dimensionally regulated electromag- 
netic form factor of a massless quark W 

r M = u(p') 7/i u(p) F(Q 2 ,e) , (1.1) 

with Q 2 = (p' — p) 2 . In QCD, the form factor F(Q 2 , e) exponentiates as 161 

, Q 2 x 

F(Q 2 , e ) = exp | - \ J ^-[K(a s (f, 2 ),e) + G{r ] 2 / f , 2 ,a s ^ 2 ),e)] 1, (1.2) 
^ o J 

where, because F(Q 2 ,e) is invariant under renormalization, the functions K and G satisfy 

dK_ = _^dG 
dfj, dpi 

Infrared singularities are summarized by K, which contains at most a single pole in e = 2—d/2, 
while G, which summarizes ultraviolet behavior, has at most a single logarithm in r] 2 /fx 2 . K 
and G are each power series in a s , which may be read off from low-order calculations.' 61 For a 
fixed coupling, \nF(Q 2 ) has double logarithms of Q 2 | 41 and double poles in e, generated from 
the explicit integrals. 

Sudakov exponentiation is much more general than this example, however. It applies as 
well to many cross sections involving a large momentum transfer, in special limits of their 
phase spaces. As we shall see, logarithms in these cross sections exponentiate after a Mellin 
or Fourier transform. We shall refer to this range of possibilities collectively as Sudakov 
resummation. 
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The connection between factorization and resummation is already illustrated by perturba- 
tive renormalization, in which the general relation of unrenormalized and renormalized Green 
functions of fields 4>i carrying momenta pi is 



G un ( Pi: M,g ) =Y[zl%/M,g(fi)) G ien ( Ph pi, g(pi)) . (1.4) 

i 

M is an ultraviolet cutoff, and g{p) and go are the renormalized and bare couplings respec- 
tively. The independence of G un from /i and G ren from M may be used to derive renormal- 
ization group equations, 

/^^ = -£7<(<M), (1-5) 

i 

in which the anomalous dimensions 7$ = (l/2)(/id/dii) InZj appear as constants in the sepa- 
ration of variables, free of explicit dependence on either \i or M. 

Very similar observations apply to structure functions in deeply inelastic scattering (DIS), 
which we may write in a convolution form, neglecting particle labels, as 

1 

F(x, Q 2 ) = [ *»- C(x/y, Q 2 /f, 2 , gfa)) /i 2 ) , (1.6) 

J y 

where C is a short-distance coefficient function, and is an infrared-sensitive parton distri- 
bution. The parameter \i here plays the role of a "factorization" scale, separating long from 
short distances. The expression (1.6) factorizes into a product under moments Jq dxx N_1 , 

F(N, Q 2 ) = C(N, Q 2 //i 2 , g(jj,)) $(N, giti)) . (1.7) 

Since F is physical and thus independent of //, we again find a simple renormalization group 
equation for C and 

^ ~ 7jv ^0 in ^ = ° = + ^^)) in ^ ' ( ls ) 

where the anomalous dimension depends on N and g, the only dimensionless parameters 
shared by the two functions on the right-hand side of (1.7). 
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We will show how Sudakov behavior can be derived from a generalization of the single- 
scale, two-function factorization illustrated by DIS above, to a two-scale, multi-function fac- 

[7J8J9] 

torization. Such a factorization is characteristic of physical quantities that describe hard 
scatterings, and which are sensitive to three generic regions of momentum space, associated 
with three sorts of quanta. These are: off-shell ("hard") partons at short distances, fast- 
moving on-shell ( "collinear" ) partons near the light cone, and long- wavelength, soft partons. 
The separation of a cross section into separate functions for each type of excitation requires 
multiple factorizations. The independence of the final expression from the details of these 
factorizations, however, leads to a consistency equation, and consequently to Sudakov resum- 
mation. This process is related to the use of matching conditions in effective field theories.' 101 

A number of phenomenologically useful and theoretically important resummations of this 
type have been derived or postulated for quantities of experimental and theoretical interest, 
including the Sudakov form factor itself,' " 1131 the transverse momentum distribution of back- 
to-back particles in e + e~ annihilation 121 1111 and of Drell-Yan pairs, 1121 the total Drell-Yan cross 
section' 131 [14 ' 1151 [161 1171 and a large class of infrared safe event shapes in e + e~ annihilation) 181 as 
well as elastic scattering amplitudes at very high energy' 191 Closely related expressions apply 
to the momentum dependence of vacuum expectations of products of Wilson lines (ordered 
exponentials) 1201 . Although all related, the derivations of Sudakov resummations for these 
quantities have been sufficiently varied in presentation and emphasis that a single derivation 
from an underlying principle should be useful. In this paper, we shall provide such a derivation 
based on a common factorization into functions sensitive to hard, collinear and soft partons. 

Yet another motivation for studying resummations of this form is that they naturally 
result in exponentials of integrals over momentum scales of the running coupling. The sin- 
gularity of the perturbative a s (fj?) at jx = Aqcd is sometimes referred to as an infrared 

[21] 

renormalon , whose presence and nature may shed light on power-suppressed corrections 
to cross sections of interest. Although we shall not investigate renormalons directly in this 
paper, the generality of Sudakov resummation should help emphasize their relevance. 

The factorization into short-distance, jet and soft functions is quite general, applicable to 
many weighted cross sections with large momentum transfers, and to some elastic amplitudes 
as well, when considered in impact parameter form.' 191 In the following section we shall 
use this factorization to derive the general Sudakov resummation formula for electroweak- 
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induced processes (Drell-Yan, deeply inelastic scattering, e + e~). This is the central result of 
this paper, from which we draw applications in Sections 3-5. Section 3 briefly discusses the 
extension of the formalism to QCD hard scattering [19] [22] , which shows that exponentiation 
applies as well to cross sections involving jets and heavy quark production. In Section 4 we 
discuss the reexpression of our result as an evolution equation in terms of momentum transfer, 
whose solution is useful for the computation of dimensionally-regularized cross sections and 
amplitudes. Finally, in Section 5 we use the Drell-Yan and deeply inelastic scattering cross 
sections as examples. 

2. Resummation from Factorization 

This section begins with a review of the basic properties held in common by all the 
cross sections (and, with slight modifications, amplitudes) to which Sudakov resummation 
applies. These properties are summarized in a rather specific but relatively simple factorized 

[7] 

expression satisfied by each cross section in a particular limit of its final-state phase space . 
This underlying factorization is illustrated in Fig 1 for e + e~ annihilation (a), the Drell-Yan 
cross section (b) and DIS (c). The figure is to be thought of as a reduced diagram, showing 
on-shell lines of the amplitude and its complex conjugate, in cut diagram notation. All lines 
are massless. 

In the figure, electroweak vector bosons of large spacelike or timelike mometum q, q 2 = 
±Q 2 , are linked to partons through two hard-scattering functions H and H*, in the amplitude 
and the complex conjugate respectively. H and H* depend only on quanta off-shell by order 
Q 2 . On-shell particles with momenta of order Q fall into two jets, J\ and J2 of collinear 
particles. We assume that the final state, represented by the cut of the diagram, approaches 
an "edge" of phase space, where all finite-energy particles are concentrated within the two 
jets. We shall also refer to this as an "elastic" limit, although the final states include arbitrary 
numbers of particles. In this terminology, the diagrams of Fig. 1 represent the most general 
configurations that contribute at leading power in Q 2 in the elastic limit.' 71 

As shown in the figure, jets may be produced at the hard scattering, as in e + e~ annihila- 
tion (Fig. la), or they may originate with an incoming parton, and hence be identified with 
a parton distribution, as in the Drell-Yan process (Fig. lb). In deeply inelastic scattering 
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(Fig. lc), we have one parton distribution and one outgoing jet. In addition, each such cross 
section involves the emission of soft particles, represented by function S. The elastic limit 
may be isolated quantitatively by introducing an appropriate weight. 

By a weight, we refer to a set of functions w n ({ki}) of the momenta of the particles in 
each n-particle final state, {ki}, i = l...n. The weighted cross section is a sum over all 
final states, each integrated over its phase space with weight w n ({ki}). For definiteness, we 
assume that the w n are dimensionless, and that they vanish in the elastic limit identified 
above. We also assume the choice of weight is consistent with the cancellation of infrared 
divergences due to final-state interactions. We shall refer to such a weight as being infrared 
safe. This requirement is met if the w n ({ki}) are symmetric in the momenta ki and satisfy 
the relations 1231 

w n (ki, ... (1 - a)k n -i, ak n -i) = w n -i(ki, . . . fc n _i) , (2.1) 

which ensures that two states that differ by the emission of zero-momentum particles, or by 
the rearrangement of momenta among collinear particles, have the same weight. 

In nth-order perturbation theory, we encounter logarithmic enhancements in the elastic 
limit on a diagram-by-diagram basis, with a™ [(ln 2n_1 w) / w] typically the most singular be- 
havior at order a™. These are the singularities we shall organize. In dimensionally- regularized 
perturbation theory, these logarithms appear from the expansion of kinematic combinations 
such as 

] ^nr, (2-2) 

with a = 1 or 2 (and e = 2 — d/2), examples of which we shall encounter below. 

Not every infrared safe weight leads to exponentiation, however . We will need to assume 
in addition that near the elastic limit, the contributions to the weight of particles within the 
jets and soft function are independent and additive, 

w = w\ + W2 + w s , (2.3) 

with corrections that vanish as w 2 for small w. 
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Consider the processes of Fig 1. For the e + e~ annihilation cross section, we may identify 
w = 1— T, with T the thrust 1181 , although similar enhancements may occur in other variables as 
well. For the inclusive Drell-Yan cross section' 131 , enhancements occur in the variable w = 1 — 
Q 2 /s = 1— r, while for incusive DIS they occur in 1— x, with x the Bjorken scaling variable' 131 . 
Finally, for the Drell-Yan cross section at measured transverse momentum [12] da /dQ 2 d 2 QT, we 
may take w ~ Qt/Q- 

In the elastic limit, we associate with the jets momenta p\ and p2, with Q 2 = ±(pi ±P2) 2 
the scale of the hard scattering. For two outgoing jets, as in e + e~ annihilation, we may 
regard each p^ as the lightlike jet momentum in the elastic limit, and Q 2 = (pi +P2) 2 - For 
an incoming jet, as in DIS, p\ may be taken as the lightlike momentum of the initial-state 
parton. In this case, Q 2 = —q 2 = —(pi — pi) 2 - 

The general factorization form near w = 0, applicable to all these cross sections is 

[ dw\ dw 2 dw s 

a(w) =H(pi/n,p2/fiXi) / 

J Wi w 2 w s 

x j x { Pl ■ Ci//i, ^i(Q//i) ai ) J 2 (P2 ■ C 2 //WQ//i) a2 ) ( 2 - 4 ) 

x S(w 8 (Q/fi),Vi,Ci) S(w -wi-W2- Wg) , 
with corrections that vanish as powers of w. The specific kinematics in (2.4) is defined by 

Pi = Q + v u p 2 = Q~v 2 , q 2 = ±Q 2 = ±2Q+Q- , (2.5) 

with Vi dimensionless lightlike vectors. As indicated, we shall take v% = 5^+, v% = 5^- for 
definiteness. The factorization eq. (2.4) often applies even if the cross section is not infrared 

[7] 

safe, but may be factorized according to the standard procedure. In Drell-Yan, for instance, 
collinear divergences are factored as usual into the jets Jj, which play the role of parton 
distributions, as noted above. In (2.4), Jj, S and H represent the contributions of the two 
jets, of the soft quanta not part of the jets, and of the hard quanta, respectively. We have 
explicitly exhibited overall factors of 1/iUj, l/w s , which we readily derive from power counting 
for the soft and jet functions in (2.4) ,' 131 

We may think of the vectors (j and £2 as gauge-fixing vectors ^-^4 = 0, used to define 
the jets Ji in terms of particular matrix elements.' 131 In such axial gauges, the factorization of 
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Sudakov double logarithms into the jets is automatic J . Although the jet definitions may be 
gauge-dependent, the cross section itself, of course, remains gauge invariant. 

The arguments Wi(Q/fj) ai represent the combination of weight and momentum transfer 
Q through which the jets depend upon the weight. The value of ai in a given cross section 
depends on the underlying kinematics. In the analysis of the x — > 1 limit of DIS in Ref. [13] 
(where w — 1 — x), for instance, a = 1 for the initial-state jet and a = 2 for the final-state jet. 
We assume (as is the case for all the examples we know of) that the soft function depends 
only on the variable w s Q. In principle, other combinations of Q and w-dependence in eq. 
(2.4) are possible, but we shall not describe them here. Because of their (^-dependence, the J{ 
depend on pi ■ Q and, in general, ( 2 . The variables pi • Q, which depend on Q, are independent 
of WiQ a \ which is defined by kinematics. 

The convolution (2.4) factorizes into a simple product under the transformation 



oo 

—Nw 



a(N) = J dwe^ w = H(p 1 /fi,p 2 / f i,Q)S(Q/fiN,v i ,QMw) 
o 

x h{pi ■ CiM Q/W 1/ai ) Mp2 ■ C 2 //i, Q//iiv 1/a2 ) , 



(2.6) 



where, for instance, J = / °° (dw/w) exp[— Nw] J. Singular behavior in the w — > limit 
is reflected by growth for large values of the conjugate variable N. This is the significance 
of the lower limit on the w integral. For large N, contributions from the upper limit are 
suppressed, and it may be extended to infinity. We consider iV as a complex variable, so 
that (2.6) represents both Laplace and Fourier transforms. For large \N\, it serves as well 
as a representation of the Mellin transform, since for w — > 0, exp [-Nw] ~ (1 - w) N , with 
corrections that vanish as w, which is the accuracy of the basic factorization formula (2.4). 

In general, the individual functions H, J{ and S each require renormalization, which we 
assume to be multiplicative. Their renormalization scale dependence is determined by a set 
of anomalous dimensions, 

/^hi# = -7fr(a a ) 

^ lnJ i = _ 7Ji(«a) ( 2 - 7 ) 
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We assume for this discussion that these anomalous dimensions are independent of the moment 
variable iV and of the £'s. The independence of cr(N) from the renormalization scale \x leads 
to 

!H + ls + ^2lJ t = 0, (2.8) 

i 

that is, the renormalization-dependence of the individual functions must cancel in their prod- 
uct. 

In summary, we are concerned here with weighted cross sections that satisfy the short 
distance-light cone-long distance convolution of eq. (2.4) in terms of functions that are renor- 
malizable according to eq. (2.7). It seems likely that (2.4) has an interpretation in the language 
of effective field theories) 101 applied to massless particles in Minkowski space. 

We emphasize that in this paper we shall not prove eq. (2.4) and its related renormalization 
conditions (2.7) for any process, nor do we wish to minimize the challenge involved in carrying 
out such a proof in any particular case. For the status of such proofs, we refer the reader to 
the literature cited above. We shall simply assume the result, and show that it is restrictive 
enough to imply very specific resummation formulas for the small-w limit. 

In the arguments below, we shall treat the choices of Q as analogous to an ambiguous 
choice of renormalization scheme. Like a choice of factorization scale, this very ambiguity 
results in strong restrictions on the cross section. 

The general factorization eq. (2.4), then, cannot depend on the precise choice of (i or (2- 
A small change of (j, for instance, leads to a relation between the first-order changes in H, 
S and Ji, which must cancel. Thus, under a variation of p\ ■ (j at fixed ( 2 (and any other 
implicit (j-dependence) in the functions, we have the relation 

= (p i < i d^~ 1 H ) Jl -^~ S + H ( w ' Cl ^5cT Jl ) j2§ + nJlJ2 ( Pl ' Cl ^5ci^) ' 

'(2.9) 

( J2 is independent of (j) or equivalently, 

Pi ' Ci o r ln Ji(pi -b/frQ/nN 1 /" 1 ) = -pi - Cits -lnH( Pl /fi,p 2 /fi,Ci) 

dpi ■ Ci dm ■ Ci (2 10) 

- Pi ■ Cit, In S(Q/fiN, Vi, d) . 

dpi ■ Ci 
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The logarithmic derivative of J\ with respect to p\ -(j may depend on the hard scale Q ~ pi -(j, 
or the soft scale Q/N. From (2.10), this dependence is additive in terms of one function that 
depends on the common variables of J\ and H and another that depends on the common 
variables of J\ and S. For H these are a s (fx 2 ) and p\ ■ (l/ f 1 - F° r S they are a s (fj, 2 ) and 
another variable that depends upon ai, given by 



i Q „ Q ai _Q' ai 



(2.11) 



where the final form serves to define Q' a . 
We may thus write (2.10) as 

Pi ■ Ci^— In Ji(pi • Ci//i, Q/ziiV 1 ^) = G ( P1 • Ci//i, a s (/i 2 )) + K(Q'JfiN, a s (fi 2 )) , (2.12) 
dpi ■ Qi 

in which the function G matches the variation of the hard part and K the variation of the 
soft part, 

G^-pih—^—lnH, 

V . (2.13) 
K = -pi ■ (i- —In S. 

dpi ■ Ci 

Eq. (2.12) is the basic consistency condition referred to in the introduction. 

Because J is multiplicatively renormalized, the combination G + K is itself a renormal- 
ization invariant 121 [3) 



, d d 



[2] 

or, again separating variables, 



G( Pl -Ci//Wm 2 )) +K(Q' ai /^N,a 8 (^)) 



0, (2.14) 



pi—K(Q' ai /Npi,a s (pi 2 )) =- lK (a s (fi 2 )), 
^ (2.15) 

^djl G ( Pl ' ^/^' as ^ 2 )) = ^(a s (/^ 2 )) , 
with jk a Sudakov anomalous dimension. This relation allows us to reexpress G + K in a 
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form which relates the scales pi ■ (i and Q' /N, 



pi-Ci 



K(Q'jNv, « s (/i 2 ))) + G(pi • ClAi, « s (/i 2 )) = ~ I % M<*a<J*' 2 )) + A'(a.((pi • Ci) 2 )) , 

(2.16) 



Q' 01 /iV 



where the functions A and A' are defined by 



d 

A(a s ) = lK (a s ) + P(g)—K(l, a s ) , 



(2.17) 



and 



A' (a 8 ) = K(l,a s ) + G(l,a s 



(2.18) 



A solution to both (2.12) and the renormalization group equation for the jets in (2.7), 
which organizes both pi ■ Q- and iV-dependence is 



J(p ■ {fa, Q/fiN^M^ 2 )) = J(l, 1, a s {Q 2 /N 2/a )) exp 



J fuM>?)) 



Q/N 1 /"- 



(2.19) 



x exp 



- / f{ / jM°.(e))-A'M\*)}], 



Q/N 1 /" Q a /X a - 1 N 



where we have used the definition in (2.11) for Q' a . 

In the following, we take p-( — Q. Shifts of this value for the form p ■ ( = CQ, with C a 
constant, may be absorbed into redefinitions of A' [2] . Finally, using the renormalization group 
equations (2.7) for H and S, and the relation (2.8) between their anomalous dimensions, we 
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find 



In a(N) = DiMQ 2 )) + D 2 (a s (Q 2 /N 2 ),a) 

Q A 

- 2 / x[ / f ^M£ 2 ))-*MA 2 )) 

Q/N 1 /"- Q a /X a - 1 N 

= D^asiQ 2 )) + D 2 (a s (Q 2 /N 2 ), a) 

Q 



-2 



Q/N 1 /" 



Q/iV 



ln(g/0^(« s (a)-^(a s (e 2 )) 



Q/iV 



ln^N/Q)A(a s (e)), 



where for simplicity we have taken a\ = a 2 = a, and where 

D 1 (a s (Q 2 )) = \nH(lA,a s (Q 2 )), 
D 2 (a s (Q 2 /N 2 ),a) = \nS(l,a s (Q 2 /N 2 )) + £ In 4(1, 1, a s {Q 2 /N 2 l a )) 

Q/(N) 1/a 



~ i 



T 



7s(« s (A 2 )), 



Q/JV 



5(a s ) = -jh + A'. 



(2.20) 



(2.21) 



Note that the last term on the right-hand side of (2.20) vanishes for a = 1, in spite of its 
overall factor of l/(a — 1). Also, we note that in cross sections involving partons in the initial 
state, such as DIS and Drell-Yan, the function D 2 will in general contain collinear singularities 
(1/e in dimensional regularization). 

At one loop in a s , denoting f(a s ) = (a s /ir)f^ + . . ., we have 



\na(N) = D^ + + (a s /n) 



+ - A<» ( -) In 2 N + 25« ( -) In N 



. (2.22) 



Thus, from the one-loop cross section the functions A and B and the combination D\ + D 2 
may be simply read off. Higher orders in the functions are similarly determined to higher 
loop order. 
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The result (2.20) shows that the cross section is an exponential, whose exponent has at 
most double logarithms of the weight variable N for fixed coupling, with higher powers of 
N generated only by expansion of the running coupling. This is the basic Sudakov expo- 
nentiation, in the form derived, for example, by Collins and Soper' 21 , and in Ref. [19]. The 
iV-dependence of couplings in the function D2 is, of course, consistent with these results. 
These formulas, however, must be restricted to values of N for which a s (Q 2 /N 2 ) remains 
well-defined. To go beyond this region, at least formally, we shall use the solution to an evo- 
lution equation in momentum transfer, which we shall derive in Section 4 below. Before doing 
so, however, we observe how these results may be generalized to purely QCD hard-scattering 
processes, such as inclusive heavy-quark or jet production. 



3. Sudakov Resummation for QCD Hard Scattering 

The methods of Section 2 above apply rather directly to processes that proceed through 

[24] [22] 

color exchange, such as heavy quark production . To our knowledge, however, until 

[22] 

recently the only treatments beyond leading logarithm in the literature for QCD hard- 
scattering processes were for the independent-scattering description of hadron-hadron elastic 

[19] [25] 

scattering, and for the scattering of Wilson lines. 

In place of eq. (2.4), we now have a somewhat more general form, in which there may be 
more than two jets, and in which the hard and soft scattering functions are labelled by their 
color content, 



a(N) = j dwe- Nw a(w) 

(3.1) 
= ^i//j(pi//i,P2//x,0)^j(Q//iA r ,^,Ci) \[Ji(pi<i/^Q/»N l l ai ). 

IJ i 

The indices IJ refer to the color structure of the hard scattering. Because soft gluons decouple 
from the jet functions, the jets are diagonal in color after factorization, just as for parton 
distributions. There are two indices in the typical cross section that describes a single hard 
scattering, one from the amplitude and one from its complex conjugate. For the high-Q 2 
elastic scattering amplitudes of protons, there are as many as three color indices, one for 
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each of three quark-quark scattering amplitudes. When the scattering involves quarks, these 
color factors may be described by singlet or octet exchange . [191[25] Quark-gluon and gluon-gluon 
scattering may be labelled correspondingly. 

[22] 

For heavy quark production, when the heavy quark mass is of the order of X1X2S, with 
Xi partonic momentum fractions, there are only two jet factors, one for each of the incoming 
partons. The heavy quark propagators involve no collinear divergences (and hence no Sudakov 
logarithms). They may thus be absorbed, in eikonal approximation, into the soft tensor Sij. 
The dynamics of the heavy quarks beyond eikonal approximation may be absorbed into Hjj. 
For two-jet production near the elastic limit, there are four J{, two representing incoming 
partons, and two representing the outgoing jets. 

The analysis of Section 2 may be repeated for (3.1), leading to an analog of the resummed 
result (2.20), except that now, because the anomalous dimensions corresponding to the hard 
and soft functions are color matrices, we find ordered exponentials for nonleading logarithms. 
To be specific, the analogs of eqs. (2.7) for anomalous dimensions are 

/^(hi#)/j = (T H (a s ))u 

/^ ln % = UiM, (3-2) 
^^InS/j = {T s {a s )) u . 

The independence of da/dw from the renormalization scale \x then leads to 



i 

We shall not work out any example' 221 of this formalism further here. Clearly, however, up 
to the matrix structure of the anomalous dimensions, the arguments of the previous section 
may be repeated. In a basis in which the color anomalous dimensions are diagonal) 19120 ' there 
is a separate exponentiation for each diagonal color structure, to leading logarithm in N for 
Sij, which is next-to-leading logarithm overall. 
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4. Evolution in the Momentum Transfer 



The Sudakov exponentiation in eq. (2.20) leads in turn to an alternate form of resum- 
mation, from its dependence on the momentum transfer Q 2 . This approach is particularly 
natural for electroweak processes, in which the hard scattering is a color singlet. 

We begin by taking the derivative of eq. (2.20) with respect to \nQ to derive the Sudakov 
evolution equation 

Q 

(N, Q 2 ) =E 1 (a s (Q 2 )) + E 2 (a s (Q 2 /N 2 ), a) - 2 J | A^ 2 )) 



Q/N 1 / a 

Q/N" 



Q/N 2 /°- 

where, (suppressing arguments of the functions), 

E 1 = P^-D 1 + 2B, E 2 = (3^D 2 -2B. (4.2) 
dg dg 

The structure of this equation is enough to imply an alternate resummed expression, of the 
sort derived, for instance, in Ref. [13], in which the exponent is expressed as an inverse 
moment. 

Eq. (4.1) shows that to leading power in N, any factorizing cross section obeys an evolution 
equation of the form 

Q2 w d{N ' q2) = W{N ' q2)HN > q2) • (4 - 3) 

To be specific we may take a(z, Q 2 ) to represent, for instance, the Drell-Yan (DY) cross 
section a DY (z,Q 2 ) or a deeply inelastic scattering (DIS) structure function F DIS (z, Q 2 ). 

The inverse transform of (4.3) is a convolution, whose precise nature depends on which 
transform is appropriate to the cross section. To be specific, and to relate to the examples we 
discuss in the following section, we choose the Mellin transform, which leads to the familiar 
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evolution equation [26] , 



W(z,Q*) = I ^N- Z W(N,Q 2 ), 



(4.4) 



where in the second line we have shown the formal inverse relation between W and W. 
Working in d > 4 dimensions, we can take these quantities to be normalized such that 



a(z,Q 2 = 0) =5(1 -z) 



(4.5) 



The physical content of this condition is that for zero momentum transfer there is no radiation 
when the scattering is through an electroweak current. In more than four dimensions, these 



[6] 



conditions are explicitly realized order-by-order in perturbation theory. 

The solution to eq. (4.3) with boundary condition a(N, 0) = 1 is simply 



a(N,Q z ) = a(N,Q z /^,a(^),e) 
Q 2 



= exp 



(4.6) 



where as above d = 4 — 2e > 4, and a is defined by 



■Is - ">(/' 2 ) 



7T 



(4.7) 



In eq. (4.6), we have made explicit the perturbative dependence of W on the renormalization 
scale \i. 

We can now employ the invariance of the physical quantity W under changes in the factor- 
ization scale ii. There is considerable freedom in how to proceed. Since the LHS of eq. (4.6) is 
a renormalization group invariant, so is the exponent, and the function W (z, £ 2 / /j? , a(jj?) , e) 
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as well. From eq. (2.20), we easily conclude that the coefficient of a(fi 2 ) M in W(N) has 
powers of In N up to ln^ +1 N. We can now show that the general form of the coefficient of 
a^(/i 2 ) in the exponent is 

M 

W ( M \z,Q 2 /^,e) = ^c ( j M) (Q 2 /i2 2 ,e)V j (z) + f (M \Q 2 //j 2 ,e)5(l - z) 

7^0 ( 4 -8) 

+ h( M \z,Q 2 /» 2 ,e), 

where h^ M \z) is regular at z — 1, and where the plus distributions T>j(z) = [ln J '(l — z)/l — z] + 
are defined as usual by 

1 1 X 

J dzg(z)[f(z)]+ = J dz(g(z) - g(l))f(z) - g{\) j dzf(z) , (4.9) 

xx 

with g(z) an arbitrary but smooth function. To show (4.8), we recall that for large N, 

I 3+1 

/ dzz N ~ 1 V j {z) = £t jfc ln fc iV, (4.10) 
{ k=o 

with corrections that fall off as 1/N, where the matrix [t] is invertible (because it is triangular; 
see for example Table 1 in 
in eq. (4.8) have moments 



see for example Table 1 in [14] for explicit coefficients). Then the plus distributions, W^f^ 1 



(4.11) 



M+l M M 

wi, M \N,Q 2 /v 2 ,a^ 2 ),e)=J2^N £ ^ + £4% 

k=l j=k-l j=0 

M+l 

-EC 1 "'". 

where = ^2cj M ^tjj~ is a function ofQ/n and e in general. 

The functions f( M ^ and h/- M ^ in eq. (4.8) are Mellin transforms of constants and terms 
behaving as powers of 1/N respectively. Since the renormalization group does not affect the 
N dependence, the fr^'s are RG invariant, and by invertibility of the matrix [t], so are the cfs. 
Then the series f(z) = XIm f^ M ^ and h(z) = XIm h^ M ^ are also RG invariant. Later we will 
use this freedom in treating parts of the resummed partonic quantities differently under the 
renormalization group. 
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For now we apply renormalization scale invariance to the full W in eq. (4.6): 



W(z, e/^ 2 , «(/i 2 ), e) = W(z, «K V, «(// 2 ), e), e) , (4.12) 

where we are free to choose v to be a function of z (and we shall below). To include all 
corrections of the form (alnN) n in the exponent, we need only up to the two-loop running 
coupling' 141 [2?1 . Here we need the dimensionally continued version of a s , which we have denoted 
a(X, a(fi 2 ), e). 

The defining equation of the ci-dimensional running a is 

xl _ e d[\*a(\,atf),e)] = _^ 2(A) a{ ^ e) _ ha3{K ^ e) (4 13) 

with boundary condition a(X = 1, «(/i 2 ), e) = a(fi 2 ). Here 62 = (HCU — 2nj)/12 and 
63 = 34(7^/48 - {20C A /3 + 4C F )n f /32. The solution, linearized in 63, is: 

x ' a ^ Q '(" 2 )- f » = i-^avU') + I ti-ffw/ '^' 1 ' (4 ' 14) 

with 7 (A«, £) = b -f (A- - 1), /(A e , a, 6) = 1 - A- " (l + tfe) ln(l - T<A e , e)a). 

Given eq. (4.12), the exponent in eq. (4.6) can now be expanded as a power series in 
a(X, a(/i 2 ), e), with A = v£ 2 /fi 2 , 

1 vQ 2 /^ 2 

a(N,Q 2 /fi 2 ,a(fi 2 ),e)=exp\ [ dzz N - 1 [ ^{a(\,a(fi 2 ),e)W^(z,l/u,e) 

H { X (4-15) 

+ a 2 (\,a(fi 2 ),e)W (2 \z,l/u,e)} 

[13] 

In this resummed expression, the exponent is itself in the form of a moment. A similar 
analysis may be given for Fourier transforms! 281291 The <i-dimensional running coupling a in 
eq. (4.15) has a usual perturbative ( "renormalon" [21] ) singularity at A = exp[— l/afe], as seen 
from eq. (4.14). Eq. (4.15) must therefore be considered as a resummation order- by-order in 
a perturbative expansion of its exponent, without further information on how to treat the 
running coupling at very low scales. 
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The functions can be determined by choosing v — 1 in eq. (4.15), and ex- 

panding it to two loops in a(/i 2 ), 



Q 2 ^2 In a(N,Q 2 /^ 2 ,a(^),e) =a(Q 2 /fi 2 , a(/i 2 ), e)W^(N, 1, e) 



+ a 2 (Q 2 /fi 2 ,a^ 2 ),e)W^ 2 \N,l,e), 



in terms of a s (fj, 2 ), using (from (4.14)) 



(4.16) 



X e a(X, a(fi z ), e) = a(fi z ) + 7(A e , e)a 2 (/i 2 ) . 



(4.17) 



We then find for the perturbative expansion of the moments 



d 



wW(N,l,e) = (Q 2 /^rQ 2 ^^\N,Q 2 /^ 2 ,e) , 



and 



W( 2 \N, l,e) = (Q 2 //^ \ ^ ^ 2 \N,Q 2 /^ 2 ,e) - \ 



a {1 \N,Q 2 /n 2 ,e) 



d 



dlnQ^ 



■<rW(N,<?/v?,e 



2 i. .2 



(4.18) 



(4.19) 



These results may be inverted to derive W^ 1,2 \z, l,e). The full functions W^ 1,2 \z,C, 2 / fi 2 ,e) 
may then be constructed by reexpanding a(£ 2 //i 2 , a(/i 2 ), e) in a(/i 2 ) in eq. (4.12), using (4.14) 
with A = £ 2 //i 2 , 



2\ 2e 



^ 2 )(,,e 2 /^e)=(^j ^(.,M)+€V)V)(|] ^W(,,l,e). 



,2\ e 



(4.20) 



According to eqs. (4.15) and (4.18) the complete one loop result exponentiates in moment 

[13] 

space , up to possible corrections that behave as 1/N. By comparison of eq. (4.15) with 
the general resummed expression in moment space, eq. (2.20), we observe that we are able 
to choose v in such a way as to absorb all logarithmic 1 — z dependence into the running 
coupling. When this is done, leading and next-to-leading logarithms in N will be resummed 
by the terms in (4.20). We shall make use of this result in the following section. 
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5. Example: Resummation for Inclusive DIS and DY 



As a final topic, we illustrate the simplicity of the resummation method of the previous 
section, by applying it directly to the Drell-Yan and deeply inelastic scattering cross sections. 
We go on to verify that the resulting expressions reproduce known finite-order and resummed 
expressions for the hard-scattering function in the Drell-Yan process in DIS and MS schemes. 

Consider first inclusive DIS, 

V(q) + h(P) -> X, (5.1) 

where V is a vector boson with spacelike momentum q, and h is a hadron. The cross section 
for this process is customarily expressed in terms of structure functions F{ [i = 1,2,3). We 
focus on the case V = 7*, and the function F 2 . It can be factorized in moment space, as 
noted in the introduction, 

1 

F 2 {h) (N,Q 2 ) = f dxx N - 1 F 2 (x,Q 2 )=J2^(N,Q 2 )<j )i/h (N,Q 2 ), (5.2) 


where Q 2 = —q 2 and x = Q 2 /2P ■ q. Here the Cj are partonic coefficient functions, and the 
parton density functions for hadron h. (In the following, we denote Mellin transforms by 
their arguments, and drop the tilde of previous sections.) We will examine F^ for large N, 
where the quark contribution i = q is dominant. Henceforth we will denote the contribution 
of a single quark to F^ by F DIS , dropping the hadron label. By itself, (5.2) does not yet lead 
to Sudakov exponentiation in moment space. In particular C is still a function of N . F DIS , 
however, obeys another factorization theorem near the edge of phase space, i.e. at large N 
(eq. (3.14) in Ref. [13]), 



F D1S (N,Q 2 ) = 



2\ 2 



(5.3) 

in terms of a modified quark distribution . p\ is the momentum of the incoming hadron 
(parton, in perturbation theory), and p 2 the momentum of the scattered quark in the elastic 
limit. In (5.3), the iV-dependence of the partonic coefficient function C q and the distribution 
4> q /h has been factorized into 0', a hard scattering function H, a jet function J, corresponding 
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to the scattered quark, and a soft function V. Each of these functions is infrared safe, but 
4>', J and V absorb all of the N- dependence, and hence singularities in the elastic limit. The 
respective N dependences of (/)', J and V may be found from eqs. (4.7), (6.1) and (7.9) of [13]. 
As above, the Q may be thought of as gauge-fixing vectors used to define the jet functions. 
The factorized expression (5.3) is of the form of eq. (2.4), and following the arguments of 
Sections 2 and 4, this alone is enough to ensure that F DIS enjoys a Sudakov resummation in 
logarithms of N. 

Analogous arguments hold for the Drell-Yan process, 

h A (pi) + h B (p2)^V(q)+X, (5.4) 

where Ha and Kb are hadrons and V is a heavy vector boson (7*, Z, ...) of timelike momentum 
q. The factorization theorem for this process reads in moment space 

1 

/ dTrN ~ 1 d% {r,Q2) = H^/A^^h/B^Q^aiji^Q 2 ) 
*J 

= E I Ho* (^) I %,A (pi ■ Cl, ^)U{^)%/B (P2 ■ C2, 7^) 



+ 0(1/N) , 

(5.5) 

where Q 2 = q 2 , r = Q 2 / S and S = (p\ +P2) 2 , and the sum is over quarks and antiquarks in 
the second line. As above, the 4>i/A{ x iQ 2 ) m the first equality are standard parton densities, 
and <7y is the corresponding partonic DY coefficient function. Furthermore, Uij(N,Q 2 ) = 
Jq 1 dzz N ~ 1 <Jij(z,Q 2 ), and similarly for the 0's. We restrict ourselves to ij = qq, which is 
the dominant production channel at large N, and denote the partonic coefficient function 
for this case by cr DY . As for -Fbis, cr DY (z) is singular when z — > 1, but satisfies a further 

[13] 

factorization , given in the second line of (5.5), in which all of the In iV-dependence of its 
moment is absorbed into jet functions ip an d a soft function U, precisely as in eq. (2.4) above. 
The explicit iV-dependeces of ip and U may be found from eqs. (5.7) and (7.8) of [13]. 

Finally, the MS distribution' 301 , which we denote 0ms, may also be treated in a similar 
fashion. This is because it differs from the F DIS as x — > 1 only by replacing the scattered 
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quark jet by an outgoing Wilson line . We shall observe below, however, certain differences 
in the resummation of this function, compared to cross sections. 

We are now ready to derive specific resummed expressions for the DY and DIS processes, 
as well as for the MS distribution. We will also determine the exponentiating two-loop 
coefficients in the hard-scattering function of DY in the DIS and MS factorization schemes. 

According to the discussion of Section 4 we need the one loop corrections to a DY , F DIS and 

[32] 

0Mg. They are in d = 4 — 2e (e < 0) dimensions 

a^(z O 2 ) - - asCp /WV r(l-e) 
a ™w>- 7r \ Q 2 J r(l-2e) 



F^(z,Q 2 /li 2 ,6) 



If 1 + Z 2 



a s C F / 47r/x 2 y T(l-e) 
~^T\Q 2 ~) r(l-2e) 



1 / I r : 2 

e 



3e/2 



1 - 



z e + z e (l-zy e (3-z) + 3z 



(5.6) 



and 



,«( z ,g2 )= a s C F /47r/i^ e r(l-e) 1 



2tt 



V Q 2 J r(i-2e) 



i + z 2 
i-2 



+ 



(5.7) 



Here Q 2 denotes the heavy vector boson mass squared in the case of DY, and — Q 2 the invariant 
mass of the photon probe in the DIS (and MS) case; fi 2 is the dimensional regularization scale, 
while z denotes the square of the ratio of the heavy vector boson mass to the parton center- 
of-mass energy in DY, and the partonic Bjorken scaling variable for DIS and MS. The plus 
distributions in the above expressions are usually expanded as 



1 




1 


— K 

+ 


ln(l -z) 


1 2 

+ 


ln 2 (l - z) 


(1-3)1+* 


+ 


_(!-*)_ 


^ (1-z) 


_ (1-*) _ 



(5.8) 



which we choose not to do here. Redefining fx 2 — ► fx 2 exp(— (ln(47r) — je)) eliminates the (47r) e 
in the expressions (5.6) , and absorbs the factors T(l — e)/T(l — 2e). 
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Let us exhibit the case of DIS in more detail, as an example. The one loop correction to 
the DIS partonic structure function leads via eq. (4.18) to Wms(z, m momentum space. 
We can write this kernel as 



Wg>(z,l,e) = 6(l-z)fg>(e) + ;f 




+ 



+ h ( DI i(z,e) , 



(5.9) 



where the coefficient functions /d IS , g^il, /ims are regular functions of their arguments at 
z = 1, given by 



We can now use (4.20) to determine the one-loop exponent, and substitute the result into the 
exponential resummed cross section (4.15). 

Given the limitations of the factorization theorem (5.3), we can discard terms of order 
1/N. Thus we drop hul(z, e) completely and expand gms(z,e) around z — 1, dropping 1/N 
terms in this expansion where convenient. 

According to the discussion below (4.10), the term 8(1 — z) and the plus distributions 
in (5.9) are separately RG invariant. We are therefore free to choose different scales v for 
these two terms in the general resummed expression eq. (4.15). The natural choice for the 
5(1 — z) /ois term is v — 1 in (4.15). Changes in \x generate terms &2hi(// / pi) at higher 
orders. The A integral may then be carried out explicitly for this term in (4.15). For the plus 
distribution term, however, the natural choice is v — 1 — z. Then, using £ = (1 — z) 1 / 2 /! in 
eq. (4.20), we absorb the factor (1 — z)~ e in (5.9) into the limit of the A integral in eq. (4.15), 
involving only the running coupling, as indicated at the end of Sec. 4. 




(5.10) 




(5.11) 



h ( D \i(z, e) = ~e[z e (l - zy e (3 - z) + 3z] . 



(5.12) 
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Treating the iV-independent term /W in this fashion, we obtain for the one-loop resum- 
mation of the DIS cross section, 



F DIS (iV,g 2 / /U 2 ,a( / x 2 ),e) 
l 

C F 

exp 



k 



(l + - 2 ) + | 



a(X, a(/i 2 ), e) 



x exp 



a(/i 2 )CF 



(11/4 + tt 2 /3) 



exp 







-^(Q 2 (l-^ 2 ,a s (/i 2 ),e) 



x exp 



a(/j, 2 )CF 



(11/4 + tt 2 /3) 



a(A, a(/i 2 ), e) 



o 



(5.13) 



where in the second form we have evaluated the 3e/2 term to lowest order in a s (Q 2 (l — z)). 
Note that as long as were are in d dimensions the A integrals in the exponent are well-defined 
at their lower limits. There is no particular problem with integrals dominated by their lower 
limits here, since F DIS is not by itself infrared safe. 

For the Drell-Yan cross section one obtains by essentially identical methods, 

l2 /,,2 



a nY (N,Q /ii ,a(/jL ),e) 
l 

= exp Cp J dz 




TV-l+e 



1-Z 
2\n /-, m 2 



:i + ^ 



j y a(A,a(/i 2 ),e) 



x exp [-a{fi 2 )C F {l/2 - tt 2 /3)] . 



(5.14) 



Here, the natural choice of scale for the plus distributions is v = (1 — z) 2 . 

It is not possible to treat the MS quark distribution in quite the same way, because 
it is not normalized at Q 2 = by electroweak current conservation. Indeed, it is usually 
regarded as a function of only a single scale, through the running coupling. Nevertheless, it 
satisfies an evolution equation' 261 of precisely the form (4.3), with Wm§ given by the standard 
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splitting functions. The problem we must therefore solve is how to fix a boundary condition 
for this equation. To this end, we note that for e < (d > 4), the dimensionally-continued 
running coupling vanishes at zero scale, a(0, e) = 0, order by order in its perturbative 

expansion in the coupling a(ii 2 ) evaluated at any nonzero scale Dimensionally-continued 
radiative corrections to fe(iV, a(Q 2 )) therefore vanish at Q 2 = 0, and we may take, just as 
for the electroweak cross sections above, 0Ms(iV, Q 2 = 0, e) = 0. The solution to the evolution 
equation is then fixed by this boundary condition, and we find 




In the first equality, T qq (N, a) are the (singular, diagonal) terms in the quark anomalous 
dimension matrix, which contribute IniV dependence. In the second equality, we have in- 
troduced an arbitrary scale fx 2 through A = fi' 2 /fi 2 , for which a(A, a(/i 2 ), e) = a(// 2 ). This 
brings our expression into the same form as eqs. (5.13) and (5.14). 

The expressions for F DIS , and a DY in eqs. (5.13), (5.15) and (5.14) are singular as 
e — > 0. However, this singular behavior is universal for all three quantities. Thus we can 
perform mass factorization to obtain the finite Drell-Yan hard part, which is a simple ratio 
in moment space: 

u qq (N, Q ,«(„)) = mg2//i2;a(/i2);6) > (5-16) 

where <f>(N) is the moment of a suitable quark distribution (mixing with other partons is 
down by 1/iV). Choosing for alternatively F DIS and 0ms m (5.16) yields the finite resummed 
hard part for the Drell-Yan cross section in the DIS and MS schemes. The factorization scale 
is set by the denominator, and we choose it to be Q 2 . Putting fi 2 = Q 2 , and returning to four 
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dimensions because u is infrared safe, we have finally 

1 



<f(iV,a(Q 2 ))=exp 



l 7V _ 1 , (i-*) 

J dz {^r=r) (1+ " 2) / f (_2CF)a(AQ 

^ (1- Z )2 



+ ^( ( 1 - ,)Q 2 ) 



(5.17) 



x exp 



«(g 2 )CF 



(9/2 + 4tt 2 /3) 



and 



^ s (iV,a(Q 2 ))=exp 



dz 



z N ~ l - 1 



1-z 



(l + z 2 ) I y(-2C F )«(AQ 2 ) 



(5.18) 



x exp 



-a{Q 2 )C F 



(1 - 2tt 2 /3) 



One may check that expansion to one loop reproduces the one loop hard parts in [32], up 
to terms of order 1/JV. We denote the iV-independent terms in the above expressions by 
A(a(Q 2 )).* Then both of the above results are of the form 



u^(N,a(Q 2 )) = A s (a(Q 2 ))I s (N,a(Q 2 )) , 



(5.19) 



where the iV-dependent exponential Is is 
I s (N,a(Q 2 )) = 



exp 



dz 



z"- 1 - 1 
1-z 



(1-z)* 



^ gi [a(\Q 2 )]+g 2 [a((l-zr^Q 2 )] 







(5.20) 



Here ms = 1, for the DIS and MS schemes respectively, and the functions g\, g 2 have finite 
expansions in their arguments: 

oo oo 

9l [a] = $>^{ n) , g 2 [a\ = ^ ■ (5-21) 

n=l n—1 

The expression (5.20) organizes all large logarithms (In A/" or ln(l — z)) in the perturbation 



expansion for u^g. Its expansion in a(Q 2 ) generates terms a n (Q 2 ) ln m N, n — 0, .., oo, m < 



* Note that they may be numerically quite important. 
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n + 1 in the exponent. To be specific, for a given order n, one can call terms for which 
m = n + 1 "leading logarithmic" , terms with m = n "next-to-leading" , and ones where m < n 
"subdominant" . In this terminology contributes at the leading logarithmic level, g^p and 

(2) (2) 

g\ at the next-to-leading level, while g 2 is subdominant. 

The form (5.20) was derived via different methods in refs. [13] and [14] for the DIS scheme, 
and for the MS scheme in Ref. [33]. The lowest order coefficients are obvious from (5.17) 
and (5.18): 

^ 1} | DIS = 2C F , g^\ Dia = -3CV/2 (5.22) 

g[ 1) \ m = 2C F , ^ 1} |ms = , (5.23) 

and agree with those found in [13], [14] and [33]. Although, strictly speaking, we have not 
proved here that the form (5.19) holds for all gf\ functions of different arguments for the 
running coupling can always be brought to the above form, by the method discussed in Ref. 
[15]. In particular, the argument of a s in #2 in eq. (5.20) is arbitrary in the MS scheme, up 
to subdominant terms, because g^lms = 0. 

(2) 

We may improve the accuracy of our resummation by determining the coefficients in 
both schemes. We do this by expanding (5.20) to order a 2 and comparing with the explicit 
two- loop results in [34], which is most conveniently done in moment space. Using methods 
described in [35], we find, as in [14], [35], 

^ 2) | DIS = = Ik 1 = C A C F (^ - C(2)) +n f C F (-f), (5.24) 

the two loop anomalous dimension first identified by Kodaira and Trentadue 1361 [3?1 W 1201 , and 



Q^U = C\ (- A + |c(2) - 3C(3)) + C A C F (-g - A C(2 ) + l^z) 

+ nf C F Q + l - C(2)) (5.25) 

9 { ?Ws = C A C F - yC(2) - ^C(3)) + n f C F (-±1 + |c(2)) . 

The numerical impact of the latter coefficients is quite small near threshold, as they contribute 
only at the sub-dominant logarithmic level. 
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6. Conclusion 



We have presented a unified method for resumming Sudakov corrections to partonic cross 
sections. We assumed the existence of a factorized expression (2.4), which is applicable to a 
large class of cross sections in the "elastic limit" of phase space, where their QCD corrections 
due to soft gluons can be large. To factorize the cross section in these regions, it is in general 
necessary to introduce dependence on auxiliary vectors Q. The freedom in choosing the Q 
leads to consistency equations which resum the large corrections. 

To illustrate the method we have rederived the large x resummed Drell-Yan hard parts 
in DIS and MS factorization schemes. These expressions are relevant to both perturbative 
resummation and power corrections. For example, in [16], the z integral in (5.19) was treated 

[38] 

by a principal value presciption for the Drell-Yan process in DIS scheme, to define the inte- 
gral in the presence of its renormalon pole [21] . It was observed that in resummed perturbation 
theory the pole manifests itself only at the level of power-suppressed ( "higher twist" ) effects. 
The precise nature of this dependence is an active subject of inquiry. Comparisons with data 
using the resummed (differential) DY cross section have also been carried out.' 171 Other ques- 
tions of theoretical interest may include the application of effective field theory methods to 
Sudakov factorization. Phenomenologically interesting applications are possible to any QCD 
reaction for which soft gluon resummation is relevant, such as high pt jet production and 
high mass heavy quark production.' 241221 We hope that the streamlined discussion of Sudakov 
resummation given above will facilitate future inquiry in this area. 
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FIGURE CAPTION 

Figure 1. (a) Cut diagram representing the factorization of eq. (2.4) in e + e~ annihilation 
in the elastic limit; (b) same for DIS; (c) same for Drell-Yan. 
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